o 

Q 



Version date: August 8, 1997 (final version) 

Ground-state magnetization curve of 
a generalized spin- 1/2 ladder 

Takashi Tonegawa-^, Takeshi Nishida-*^ and Makoto Kaburagi^ 

^Department of Physics, Faculty of Science, Kobe University, Rokkodai, 

Kobe 657, Japan 

<D '^Department of Informatics, Faculty of Cross- Cultural Studies, Kobe Uni- 
^ ■ versity, Tsurukabuto, Kobe 657, Japan 

Employing a method of exact diagonalization for finite-size systems, 
we investigate the magnetization curve in the ground state of an antifer- 
romagnetic spin-1/2 ladder with additional exchange interactions on di- 
agonal bonds, which is equivalent to an antiferromagnetic spin-1/2 chain 
with bond-alternating nearest-neighbor and uniform next-nearest-neighbor 
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^ interactions. It is found that a half-plateau appears in the magnetization 

^H ■ curve in a certain range of the interaction constants. This result is dis- 

^ cussed in connection with the necessary condition for the appearance of 
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the plateau, recently given by Oshikawa et al. 

Key Words: ground-state magnetization curve, generalized spin-1/2 lad- 
der, spin-1/2 chain with competing interactions, half-plateau 

Address for Further Correspondence: 
Takashi Tonegawa 

Department of Physics, Faculty of Science, Kobe University, 

Rokkodai, Kobe 657, JAPAN 

Tel.: +81-78-803-0541 Fax: +81-78-803-0722 

e-mail: tonegawa@kobe-u.ac.jp 



There has been a considerable current interest in the study of quantum 
spin systems with competing interactions, which exhibit a variety of fasci- 
nating phenomena originating from frustration and quantum fluctuation. 
In this paper we investigate the magnetization curve in the ground state of 
an antiferromagnetic spin- 1/2 ladder with additional exchange interactions 
on diagonal bonds. We express the Hamiltonian describing this system in 
an external magnetic field as 
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where Sj^i is the spin-1/2 operator at the jth leg and the ith rung; l+a 
and 1 — a (0.0<q;<1.0) are, respectively, the interaction constants between 
neighboring spins along the rung and diagonal bonds; J {J > 0.0) is the 
interaction constant between neighboring spins along the leg bond; H 
{H > 0.0) is the magnitude, in an appropriate unit, of the magnetic field 
applied along the z-axis; N is the total number of spins in the system and is 
assumed to be a multiple of four. We impose periodic boundary conditions 
{Sj^]S[^i = Sj^i). It should be noted that this system is equivalent to an anti- 
ferromagnetic spin-1/2 chain with both bond-alternating nearest-neighbor 
interactions, the interaction constants being l+a and 1 — a, and uniform 
next-nearest-neighbor interactions, the interaction constant being J, which 
compete with each other. 



The ground-state [1] and thermodynamic [2-4] properties of the sys- 
tem in the case of a = 0.0 and H = 0.0 have been studied extensively. In 
particular, it has been found that as the value of J increases, the phase 
transition from the massless spin-fluid phase to the massive dimer phase 
occurs at J =Jc~ 0.2412 in the ground state [4-6]. As is well known, the 
exact ground-state magnetization curve in the case of a = J = 0.0 has been 
obtained by Griffiths [7]. The ground-state magnetization curves in the 
case of a = 0.0 [8] and in the case of J = 0.0 [9] have also been numerically 
calculated for various values of J and a, respectively. 

Obtaining the ground-state magnetization curve of the present system, 
we employ Sakai and Takahashi's method [10], by the use of which they have 
discussed the ground-state magnetization curve for a uniform spin-1 chain. 
The outline of this method, which is based on a method of exact diagonal- 
ization for finite-size systems combined with the conformal field theory, can 
be summarized as follows. Let Eq{N,M) be the lowest-energy eigenvalue, 
within the subspace determined by the value M = ^- ^ X^£=i '^j li^^^ 1^ 
• • •, N/2), of the Hamiltonian Tio for a given N. Then, the conformal field 
theory predicts that, if the state with the energy Eo{N, M) is massless, the 
asymptotic behavior oi Eo(N, M) in the thermodynamic {N ^-oo) limit has 
the form [11] 

^°^^l^~e(m)-C{m)^ (N^oo) (2) 

with m = M/N, where £{m) is the lowest energy per spin for a given m in the 
thermodynamic limit and C{m) is a positive constant which is proportional 
to the product of the central charge and the sound velocity. Minimizing 
with respect to m the total energy, £tot=£{'>in) — Hm, per spin of the system 
described by the Hamiltonian TL, we can obtain the equation which relates 



the ground-state magnetization im) per spin in the thermodynamic hmit 
with H as 

e'({m)) =H. (3) 

From Eq. (2) we can derive 

AEo(iV;M,M-l)^£'(m_o)-^£"(m_o)^ (TV ^ oo) , (4a) 

A£;o(iV;M + l,M)^£'(m+o) + ^£"(m+o)^ (TV ^ oo) , (4b) 

where AEoiN; M, M -l) = EoiN, M)- EoiN, M -1). The essential point 
of Sakai and Takahashi's method [10] is that as far as the massless states are 
concerned, Eqs. (4a, b) hold when A^ is sufficiently large, and thus the value 
H = £'{m) [see Eq. (3)] of the magnetic field for a given m (0<m<l/2) 
can be estimated by making an extrapolation which uses these equations 
to be H = £'{m-o)=£'{mj^o). For the estimation of Hco^s:'{0), which is 
the critical field at which {m) starts to increase from zero in the ground- 
state magnetization curve, we apply Shanks' transformation [12] to the 
sequences {A£'o(A^; 1,0)}, following Sakai and Takahashi's procedure [13]. 
Furthermore, the saturation field H^ = £'(1/2) can be obtained analyti- 
cally, since it is straightforward to diagonalize the Hamiltonian Tio within 
the M=(7V/2)-l subspace; it is given by i/s = Max(2, 1 + 2 J+ a) in the 
present case. Recently, we have successfully applied the method to the 
case of an antiferromagnetic spin-1 chain with bond- alternating nearest- 
neighbor interactions and uniaxial single-ion-type anisotropy [14]. 

We obtain the ground-state magnetization curve in the thermodynamic 
limit for J = 0.1, 0.2, 0.3, and 0.4, for each of which various values of a 
are chosen. In the calculation we numerically diagonalize the Hamiltonian 
Tio, using the computer program package KOBEPACK/S [15], to calcu- 
late Eq{N,M) for N = 8, 12, • • •, 24. Then, we can make the analysis for 



m=l/12, 1/8, 1/6, 1/4, 1/3, 3/8, and 5/12. Our calculation shows that 
both AEo{N;M + l,M) and AEo{N;M,M -1) are almost linear func- 
tions of 1/N at least for m = 1/8 and 3/8 in accordance with the forms 
given by Eqs. (4a,b). The values of £'(m_o) and ^'(m+o) for these m's can 
thus be estimated. In a similar way we estimate £'{m-o) and £'(771+0) for 
m = 1/12, 1/6, 1/3, and 5/12, assuming Eqs. (4a, b), although only two 
data are available. All the obtained results show that £'(m_o) and £'(7714-0) 
coincide with each other within the numerical error. For m = 1/4, on 
the other hand, Eqs. (4a, b) do or do not hold depending upon the values 
of a and J. The case where Eqs. (4a, b) do not hold is the case where 
the state with m= 1/4 is massive and therefore £'(l/4_o) is smaller than 
£'(1/4+q). This means that, in the magnetization curve, there appears 
the half ((771) = 1/4 )-plateau with the critical field Hd =£'(l/4_o) at which 
the plateau starts and that i7c2 = £''(l/4+o) at which it ends. We estimate 
the former and latter critical fields by applying Shanks' transformation [12] 
to the sequences {AEo{N;N/4:,N/A-l)} and {AEo{N; N/4: + l,N/A)}, re- 
spectively. 

We find that the half-plateau appears in the magnetization curve when 
0.5 < a < 0.95 for J = 0.1, when 0.2 < a < 0.85 for J = 0.2, when 0.0 < a < 0.8 
for J = 0.3, and when 0.0 <q;< 0.75 for J = 0.4. As an example, we depict 
in Fig. 1 the magnetization curve with the half-plateau, obtained for J = 0.2 
and q; = 0.5. Plotting versus a the critical fields Hco, Hd, and Hc2 as well 
as the saturation field Hs, we can draw the ground-state phase diagram on 
the H versus a plane; the result for J = 0.2 is shown in Fig. 2. 

We also calculate numerically the eigenfunctions of the lowest- and 
second-lowest-energy states within the M = N/4: subspace for the finite- A^ 
systems, and find that at least for a set of J and a giving the plateau. 



the lowest-energy state for M = N/A in the thermodynamic hmit is doubly 
degenerate, one of the eigenfunctions of which has the periodicity n = 4 (in 
units of the lattice constant) concerning the translational symmetry. This 
result is consistent with the necessary condition n(S'— (?7i)) = integer for 
the appearance of the plateau with the magnetization {m) {S is the mag- 
nitude of spins), which has recently been given by Oshikawa, Yamanaka, 
and Affleck [16]. Finally, it should be noted that very recently Totsuka [17] 
has clarified in an excellent way the mechanism for the appearance of the 
plateau in the present system, using a bosonization technique. According 
to this work, it becomes clear that the next-nearest-neighbor interaction 
plays a crucial role in the appearance of the plateau. 
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Figure Captions 

Fig. 1. Ground-state magnetization curve in the thermodynamic hmit ob- 
tained for J = 0.2 and q; = 0.5. The closed circles show plots of (m) versus 
H = £'{{m)). The solid lines are guides to the eye. 

Fig. 2 Ground-state phase diagram on the H versus a plane in the ther- 
modynamic limit for J = 0.2. The closed circles show plots versus a of the 
critical fields i^co, Hd^ and Hc2 and also of the saturation field H^. The 
solid lines are guides to the eye. The magnetization {m) is given by {m) =0, 
1/4, and 1/2 in the regions A, B, and C, respectively. In the remaining 
region {m) increases continuously as H increases. 
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